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ABSTRACT 

We claim that the discrepancy found between theoretical predictions for the 
cosmic mass function and those found in numerical simulations is due to the 
fact that in deriving the former all mass elements are assumed to be at the 
center of the object they belong to (the all-mass-at-center problem). By an 
appropriate treatment of this problem, using the spherical collapse model (which 
is not a bad approximation in the high mass limit), we obtain both the high mass 
behaviour found in simulations and the true assymptotic behaviour of the mass 
function (for arbitrarily high masses). Therefore, we conclude that by combining 
ellipsoidal dynamics with a suitable treatment of the all-mass-at-center problem 
(that we will show in a follow up work) a theoretical prediction for the cosmic 
mass function in full agreement with simulations may be obtained. 

Subject headings: cosmology: theory — dark matter — galaxies:clusters:general 
— galaxies: halos 



1. Introduction 



Cosmic mass functions corresponding to various cosmological models and collapsed ob- 
ject defining criterias have been accurately determined (Jenkins et al. (2001), Warren et al. 
(2006)) by means of numerical simulations covering five orders of magnitude in mass. On the 
other hand, theoretical ideas concerning these functions have developed considerably from 
those contained in a seminal work by Press & Schechter (1974). In this work, the simple 
dynamics of the spherical collapse model was combined with Gaussian statistics to derive an 
expression for the cosmic mass function. However, one of the criticism that can be made to 
this derivation is that mass elements lying at the center of underdense (or not sufficiently 
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overdense) density fluctuations on a certain scale are assigned to collapsed objects smaller 
than this scale, neglecting the possibility that these elements could lie at the center of a 
sufficiently overdense fluctuation on a larger scale (so that it should be assigned to a col- 
lapsed object larger than this latter scale). To address this difficulty, currently referred to 
as the cloud-in-cloud problem. Peacock & Heavens (1990) and Bond et al. (1991), using 
for Gaussian fields ideas advanced by Epstein (1983) for Poissonian point distributions, 
developed the, so called, excursion set approach, whereby every mass element is assigned to 
an object with the mass corresponding to the largest sphere centered at this element and 
with mean inner overdensity such that accordingly with the spherical collapse model it has 
just virialized. 

Even with this last refinements, theoretical predictions for the mass function did not provide 
a perfect fit to the results found in numerical simulations (Gelb & Bcrtschinger (1994)). 
This prompted questioning the suitability of the spherical collapse model (Barrow & Silk 
(1981)). To deal with this issue Bond & Myers (1996) used ellipsoidal dynamics within a 
general framework for identifying in the initial conditions the points at which objects of dif- 
ferent masses will form, and following their motion. Sheth, Mo & Tormen (2001) discussed 
a simple procedure for incorporating ellipsoidal dynamics into the excursion set formalism, 
and developed it to fuller extent in Sheth & Tormen (2002). Their prediction for the mass 
function is in much better agreement with those found in numerical simulations, specially 
in the low mass limit where it seems to give the correct logarithmic slope. However, in the 
high mass limit, where the functional behaviour of their mass function is dominated by the 
same exponential factor that dominates Press-Schechter expression, namely: 



with a equal to 1 (cr(m) being the rms linear density fluctuation as a function of mass), 
the discrepancy is larger than a factor 10 for the largest masses available ( ~ 5 ■ lO^^h^^ Mq 
). A value of a around 0.75 in factor (1) with Sc ~ 1.7 is in much better agreement with 
numerical results. 

This discrepancy poses a rather vexatious question which is not easy to justify theoretically. 
If collapsed objects are defined as spherical overdensities with density contrast A and 5c 
is its linear counterpart through the spherical collapse model (including shell-crossing non- 
radial motions and radial velocity dispersion), Sheth, Mo & Tormen (2001) procedure leads 
inescapably to (1) with a = 1 for the leading factor of the high mass limit. This is so because 
for the largest masses the collapse is almost spherical and the expression obtained using the 
spherical collapse model must be recovered. 

Sheth & Tormen (2002) derivation of the mass function has certain amount of uncertainty; 




(1) 



-3- 



contained in the way in which they modelled the process of virialization of the proto-objects, 
where massive shell-crossing takes place. However, this model has been built so that for 
spherically collapsing patches the actual density contrast, 6, takes the value A when the 
linear density contrast , Si, takes the value Sc- So, the above result is not affected by the 
mentioned uncertainty. 

For A ~ 200 the value of 6c is around 1.7 (1.686 is the conventional value for Einstein-deSitter 
cosmology). In this mentioned above, by comparing theoretical predictions with 

simulations one obtains a ~ 0.75. However, even if we were to use the standard spherical 
model (neglecting the dynamical ralentization due to shell-crossing and other factors), 5c 
must be close to 1.6, so that a must be around 0.85 to agree with simulations. Now, any 
refinement in the dynamical model used (i.e. including substructure) leads to a larger value 
of 5c)therefore, the fact that a is close to 0.75, cannot be explained in this manner. So, we 
must search for its explanation in another assumption of the formalism. 
In this work we show that the failure of theoretical predictions to render the correct value of 
a is due to the fact that in the excursion set formalism (as well as in that of Press-Schechter) 
the criteria for deciding to objects of which mass must be assigned a mass element is based 
upon the global behaviour of patches (spherical or eUipsoidal) centered at that element (all- 
mass-at-center) . Obviously, of all mass elements making up an object only one is exactly 
at the center, but assuming for simplicity every element to be at the center of the objects 
it belongs to is not, generally, a bad approximation. However, as we shall show, for large 
masses it leads to an underestimation of the mass function. 

We show here that if, instead of assuming every mass clement to be at the center of its object 
(with mass m), we merely assume that it is closer than Q{m) {Q being the Lagrangian radius 
of the object) to the real central element, the value of a ~ 0.75 comes out nicely. 
It must be noted that in the algorithm of Bond & Myers (1996) for obtaining mass functions 
from realizations of initial conditions, the problem is implicitly treated like here, but not 
explicit analytical treatment has so far been available. Furthermore, the question of all- 
mass-at-center has not been investigated in regard with the present problem. If Bond & 
Myers algorithm were used to this purpose, the results shown here would also be found. 
In this work we only consider spherical dynamics, since the effect considered is only important 
for large masses where the collapse is almost spherical. This simpler dynamics helps making 
more patent the essence of the effect, which shall be treated to full extent using triaxial 
collapse in a follow up paper. 
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2. Dealing with the all-mass-at-center problem 



Let F{m) be the fraction of mass making up objects with mass larger than m, and n(m) 
the number density of objects within a unit mass interval around m. The following identity 
must obviously hold: 



This relationship, in various forms (using Ln(m), o"(m), Ln(cr(m)) as independent vari- 
able), is used in most theoretical derivations of the mass function, in particular, in the 
excursion set and Press-Schechter formalisms. It is the approximations used to obtain F(m) 
that make the difference between various approaches. For example, the dynamics can be 
modelled by the spherical collapse model or by the more refined ellipsoidal model. Here we 
adopt the spherical model since the effect we intend to explain, namely, the "anomalous" 
large mass behaviour of the mass function does not depend significantly on the dynamical 
model used. 

In the excursion set approach with spherical dynamics F(m) is equated to the probability 
that for a randomly chosen mass element the largest sphere centered at this element with 
an inner linear density fluctuation equal to 5c has a radius larger than Q{m). That is, as 
mentioned in the introduction, the computation of n(m) proceeds as if all mass elements 
were at the center of the object they make up. This assumption is common to all approaches 
using expression (2) developed so far. The peak approach (Peacock & Heavens (1990)) does 
not make this assumption but it has other problems that we shall discuss in a follow up 
work. 

The aim of this work is to show that an approach using expression (2) may easily be devel- 
oped where F(m) is computed applying the spherical collapse model only to spheres centered 
at the central point of each object, and then show that the resulting mass function has the 
correct large mass behaviour. 

If the spherical collapse model applies, a randomly chosen mass element belonging to an 
object of mass m must lie, in the initial conditions, within a distance (lagrangian) Q(m) 
from a point such that a sphere of radius Q(m) centered at it has an enclosed mean linear 
density contrast, 5i, equal to 5c. 

Let us now obtain the probability distribution for 5i within a sphere of radius Q(m) centered, 
not at the center of the proto-objcct but at a randomly chosen point within it. When such 
point is at a distance (Lagrangian) q from the center, the Gaussian bivariate distribution for 
the hnear fields in scale Q(m) at q and at g = (the center) leads straightforwardly (Patiri 
et al. (2006), expression (24)) to the following expression for the distribution of 5i at q, 
— ^i{q) — ) conditioned io 5i{q — — 5c. 
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where | (Jfc P stands for the hnear power spectra. S{q,Q) may accurately be approxi- 
mated by: 

S{q, Q) ~ e-^(^)'; c = c(m) = -Ln{S{Q, Q)) (4) 

To obtain the conditional probability distribution ,Pi{5i\5i{q = 0) = 8c,Q), for 5i (on 
scale Q(m)) at a randomly chosen point within the proto-object given that 5i{q = 0) = 5^ 
we multiply expression (3) by the probability distribution for q at a randomly chosen point 
within Q and integrate over q: 



Pi{5i\5M = 0) = 5,,Q) = P{5i{q)\5i{q = 0) = S,,Q)-^dq (5) 

When 5i{q — 0), instead of taking the value 5c takes a larger value, the resulting collapsed 

object (centered at g = 0) would have a mass, m/, somewhat larger than m. So, the 
probability distribution, P2{5i\m/ > m), for 5i (on scale Q{m)) for a point chosen at random 
within a proto-object with mass ml larger than m (that is, for a point lying at a distance 
smaller than Q{ml) from the center of their proto-objects), may be obtained by multiplying 
the conditional distribution given in expression (5) (with 5i{q = 0) as variable in the place 
of 6c in (3)) by the probability distribution for di{q = 0), integrating for di{q = 0) > Sc, and 
dividing by the probability that 5i{q = 0) > Sc- 
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It is clear that expression (6) cannot be exact, because in expression (5) mf — m is 
assumed, but for Si{q = 0) larger than 5c mf is somewhat larger than m. However using 
(6) leads only to a minute error that shall be discussed in a future work. Carrying out the 
integral over Si{q — 0) in (6), we find: 

P2{0i\mf > m) = . r (7) 



27ra(m)ler/c(^) 



where use has been made in expression (3) of the approximation given in (4) and where 
the variable q has been changed to u (= ^). 

When 6i is larger than 6c for a sphere of radius Q{m) centered at a randomly chosen mass 
element, it is clear that this element must belong to an object of mass larger than m. So, 
the probability distribution for 5i within a sphere of radius Q{m) centered at a randomly 
chosen mass element must be, for 5i > 5c, equal to the probability distribution for 5i for 
elements chosen at random amongst those belonging to objects with mass larger than m (i.e. 
P2{5i\m/ > m)) multiplied by the probability that a randomly chosen element belongs to an 
object larger than m. But this later probability is simply the fraction of mass belonging to 
objects with mass larger than m, which we represent by F{m). We may then write: 

exp ( — 



F{m)P2{5i\mi>m)^ - _ Mrn^ ^^^^^ 



Using expresion (7) we find: 



From expression (9) we sec that is simply the ratio between the probability that 
5i > 5c within a sphere of radius Q{m) centered at a randomly chosen mass element and the 
larger probability that a random mass element belongs to an object of mass larger than m. 
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3. The Mciss Function 

Using in expression (2) F{m) as given by approximation (9) we find that the mass 
function, n{m), may be written in the form: 

[1 - Z(m)] 



n m) = - - — ^^exp - — ^ — ^— — ^-^ (10 
^ ^ L Vtt/ ma(m)2 2a^{m)J dml V{m) ^ ^ 



Z(m)=(-Y^^e^v( )erfc( \ dLnVim) 

In figure (1) we show c(m) (see expression (4)), V{m), Z{m) for a CDM power spectrum 
(BBKS 1986) with as = 0.9, h = 0.7,0^ = 0.3 and with 4 = 1-676. It is apparent that 
Z{m) is small compared with 1 for all masses (its largest value is about 0.26). So, expression 
(10) is essentialy the normalized Press-Schechter mass function divided by V(m). Therefore, 
it is this function that makes the most difference between the full mass function given by the 
spherical collapse model and the approximations to this function provided by either Press- 
Schechter or excursion set formahsms. In the small mass hmit V(m), (which is by definition 
smaller than one for any mass), goes to one: it is larger than 0.9 for masses smaller than 
lQi^^h~^ Mq and it is still larger than | for 3 • lO^^/i"-*^ Mq. It is only for larger masses that 
it becomes a steep function of mass, going down to 0.06 for 5 • lOt^^h"^ Mq. So, it is at the 
high mass hmit (m > 10^^/i~^ Mq) where neglecting the all-mass-at-center problem leads to 
substancial errors. 



From expression (10) and figure (1) it is apparent that the difference with the normal- 
ized P-S mass function (first factor in (10)) is considerable only for m > 10^^/i~^ M©, giving 
this expression a substantialy enhanced density for the largest masses that have collapsed at 
present. For smaller masses the difference between both mass functions is small. However, 
expression (10) must give smaller densities than P-S in the small mass limit to compensate 
for the high mass excess, since both mass functions are normalized. This is actually the case 
for masses smaller than roughly 10^^/i~^ M©, the difference increasing slowly with decreasing 
mass and stabilizing around a 20% for masses below 10^°/i~^ Mq. This difference is mostly 
due to Z{m), {Z{l{f'h~^ Mq) — 0.25), since V{m) is close to one for these masses. 

The interest of the low mass limit is mostly academic, since ellipsoidal collapse plays an 
important role in this limit. In the high mass limit, however, the collapse is close to spherical, 
so that the ratio between expression (10) and P-S (i.e. the second factor in (10)) should be 
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basically the same as that between the actual mass function and the mass function obtained 
using ellipsoidal dynamics but with the excursion set approach. To support this conjecture, 
we show in figure (2) the ratio between expression (10) and P-S (normalized) along with 
the ratio between Warren et al's (2006) fit to numerical results and Sheth & Tormen mass 
function (a = 1) (Sheth, Mo & Tormen (2001)). The basic agreement between both ratios 
strongly suggest that the appropriate treatment of the all-mass-at-center problem is all that 
one needs to add to Sheth & Tormen derivation to obtain a mass function in full agreement 
with simulations. 

The ratio between expresion (10) and P-S, for the mass range considered, is essentialy given 
by V{m)~^ (expression (9)). We see, from figure (2), that for m < 3 • 10^^h~^ Mq, V{m) 
may be approximated by a function of the form given in expression (1) with a ^ 0.25. 
So, in the high mass limit, expression (10) is dominated by a term of the form (1) with 
a ~ 0.75. The fact that a is smaller than one, that we called "anomalous assymptotic 
behaviour", is therefore explained by the approximately exponential behaviour of V(m). 
It must be noticed, however, that V{m) cannot take the form given in (1) for very large 
masses. In fact, it may easily be shown that in the arbitrarily high mass limit (as opposed 
to the hmit of the highest masses available, that we have so far considered) the dependence 
of V on a{m) is not exponential, but potential: 



So, in the truly assymptotic limit, expression (10) behaves as expression (1) with a = 1, 
which, as we mentioned in the introduction, is the value that a must take in this limit. 
Using (11) in (10) leads to the same assymptotic expression as using, in this limit, the peak 
approach, as we shall show in detail in a forthcoming paper. 

Summarizing, the high mass behaviour of the cosmic mass function observed in numerical 
simulations, although dominated by a factor of the form given in (1), do not yet shows the 
real assymptotic behaviour. The apparent assymptotic behaviour led by expression (1) with 
a ~ 0.75 results from the approximate exponential behaviour of V{m) in the relevant range 
of masses. The real assymptotic behaviour of V{m) is potential with respect to cr(m), re- 
sulting in a true assymptotic behaviour for n{m) which takes the form (1) with a = 1; as it 
must be. The function V{m), which makes the essential difference with previous approaches, 
originates from an appropriate treatment of the all-mass-at-center problem. 
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Fig. 1. — c{m) defined in (4) is given by tlie dotted line; V{m) (defined in (9) is given by the solid 
line; and Z(m) (expression (10)) is given by the dashed line. 
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Fig. 2. — Ratio between the mass function as given by expression (10) and as given by the normal- 
ized Press-Schechter expression (solid line); an approximation to this ratio given by expression (1) 
with a = 0.25 (dotted line); and the ratio between Warren et al's (2006) fit to numerical results and 
Sheth &; Tormen mass function (a = 1) (dashed line). All theoretical expressions use 5c = 1.676. 



